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ABSTRACT 


This  research  memorandum  describes  a  data  processing  technique  for 
analyzing  the  Doppler-shift  properties  of  a  radio  channel.  A  computer 
program  has  been  developed  that  computes  the  power  spectrum  of  power 
scattered  by  the  radio  channel  as  a  function  of  Doppler  shift.  This 
program  was  developed  in  conjunction  with  phase-stable  CW  experiments 
which  were  performed  on  two  HF  paths .  The  experiments  were  part  of  a 
simulation  experiment  conducted  by  Stanford  Research  Institute  for  the 
United  States  Army  Electronics  Laboratories  and  the  Defense  Atomic 
Support  Agency.1 

The  program  computes  the  power  spectrum  of  the  time  variations  of 
the  response  of  the  radio  channel  to  a  CW  tone  relative  to  the  trans¬ 
mitted  frequency. 


References  are  given  at  the  end  of  the  memorandum. 
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I  INTRODUCTION 


In  a  study  recently  conducted  at  Stanford  Research  Institute,  the 
natural  radio  channel  is  viewed  as  a  randomly  time-varying  linear  filter.2 
The  purpose  of  this  study  was  to  develop  mathematical  channel  models  that 
would  enable  one  to  analyze  the  performance  of  various  communication 
systems  in  the  presence  of  signal  perturbations  induced  by  the  propagation 
medium . 

In  this  study  the  radio  channel  is  represented  by  a  time-variant 
channel  transfer  function,  H(t,f),  which  is  a  complex  random  quantity 
that  is  wide-sense  stationary  in  both  time  and  frequency.  This  function 
is  simply  defined  as  follows.  If.  the  input  to  the  channel  is 

f  i2rr(f  4-  f  )t. 

x(t)  =  Re  •! e  j-  (la) 

then  the  output  is 

r  i2n(f  4-  f  )t. 

z(t)  =  Re  iH(t,f)e  '  j-  .  (lb) 


Complex  envelope  notation  is  employed  in  Eq .  (1);  the  symbol 
denotes  the  real  part  of  the  complex  quantity  in  the  braces, 
some  reference  frequency. 

The  second-order  statistical  variations  of  H(t,f;  are  summarized 
by  the  time-frequency  autocorrelation  function: 


Re  {} 
and  f0  is 


HH(0f,3) 


r 

=  e  : 


H*(t,f) 


H(t  4-  a,f  4- 


»] 


(2) 
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where  the  symbol  J  denotes  the  mathematical  expectation  of  the 
quantity  in  the  brackets.  The  Fourier  transform  of  R^Ca,  is  called 
the  channel  scattering  function: 


S(X,T)  =JJ  e  l2TT(X*  ‘  T0)  dadg  .  (3) 


The  channel  scattering  function  is  the  density  of  power  scattered  by 
the  channel  to  the  Doppler  frequency  X  and  the  time  delay  T.2  Thus 
the  channel  scattering  function  yields  an  immediate  physical  charac¬ 
terization  of  the  frequency-  and  time-dispersive  properties  of  the 
propagation  medium.  The  function 


T(t)  =  J'  sa,T)  dX 


(4) 


is  called  the  channel  delay  profile,  since  it  yields  the  density  of 
power  scattered  by  the  channel  as  a  function  of  time  delay.  Similarly, 
the  function 


D(X)  =  J  S(X,T)  dT 


(5) 


is  called  the  channel  Doppler  profile,  since  it  yields  the  density  of 
power  scattered  by  the  channel  as  a  function  of  Doppler  shift. 
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II  THE  CHANNEL  DOPPLER  PROFILE 


-  4-=* 


The  channel  Doppler  profile  is  the  power  spectrum  of  the  wide-sense 
stationary  time  variations  of  the  channel  transfer  function.  From  Eqs. 
(3)  and  (5)  one  can  deduce  that  the  Doppler  profile  can  be  written  as 
the  Fourier  transform  of  the  autocorrelation  function  R^Cor^O).  Observe 
that  this  autocorrelation  function  ignores  the  wide-sense  stationary 
frequency  variations  of  the  channel  transfer  function. 


(t,f)  H(t 


+  a,f)  =  RH(»,0) 


(6) 


D(X)  =  J*  RjjC^O)  e  l2nXa  dor 


(7) 


Hence  the  Doppler  profile  is  merely  the  power  spectrum  of  the  time 
variations  of  the  channel  response  (amplitude  and  phase)  to  a  CW  tone 
at  a  fixed  frequency.  It  is  evident  that  the  Doppler  profile  can  be 
estimated  with  the  amplitude  and  phase  measurements  obtained  from  a 
phase-stable  CW  experiment,  since  such  an  experiment  measures  H(t,f). 

It  is  convenient  to  introduce  the  quadrature  components  of  H(t,f): 


H(t,f )  =  U(t)  +  iV(t) 


(8) 


In  Eq .  (8)  the  frequency  dependence  of  the  quadrature  components  is 
suppressed  because  the  frequency  f  will  remain  fixed  in  all  arguments 
that  follow.  The  autocorrelation  function  R^Cc^O)  can  be  expressed  in 
terms  of  the  autocorrelation  functions  associated  with  the  quadrature 
components : 
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Y“<0>  =  R^ta)  +  »n(«)  +  i  {Ruv«y)  -  R^te)}  (9) 

where 

Ruu(a)  =  E[u(t)u(t  +  a>] 

R^a)  =  E^U(t)  V(t  +  a)] 


and  similarly  for  R  (of)  and  Ryy(o') .  In  most  cases  of  practical  in- 
2  * 

terest. 


Ruu^  Rw(a) 


Ruv(o)  =  -  RvuCo> 


which  implies  that  the  autocorrelation  function  R^(q^O)  can  be  written 


^(<*,0)  =  2Ruu(a)  +  i2Ruy(ff)  .  (10) 


By  using  Eq.  (7)  and  the  fact  that  R^Cff)  is  an  even  function  of  a  and 
RyV(o)  is  an  odd  function  of  a ,  the  Doppler  profile  is  obtained  in  the 
form 


D(\)  =  2jRuu(o')cos2n\odo' +  2jRuv(a)sin2nXckdo?  .  (11) 


$ 

These  conditions  are  necessary  for  wide-sense  stationarity  of  the 

_  -x  i2n(f  +  fn)t\  .  .  . 

process  Re|H(t,f)e  u  J-  m 
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Observe  that  the  Doppler  profile  is  symmetrical  about  zero  if  and  only 
if  RyyCa)  is  zero  for  all  a.  When 


x(t) 


Re 


{6 


i2rr(f  +  f 


0 


is  the  input  to  the  channel,  the  (two-sided)  spectrum  S  (u)  of  the  re- 

z 

suiting  narrow-band  output  process 


z(t) 


=  Re 


i2n(f  +  fQ)t 
e 


} 


can  be  expressed  in  terms  of  the  Doppler  profile 


Sz(u)  =  \  D(u  -  {fQ  +  f})  +  |  D(-  u  -  {fQ  +  f})  .  (12) 


Thus,  if  D(X)  =  0 
frequency-shifted 


for  |\|  >  f  +  fQ,  the  Doppler  profile  is  merely  a 

version  of  the  narrow-band  spectrum  S  (u) . 

z 
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Ill  DESCRIPTION  OF  THE  POWER  SPECTRUM  PROGRAM 


The  input  to  the  power  spectrum  program  that  computes  the  channel 
Doppler  profile  consists  of  two  data  sets;  each  data  set  corresponds  to 
N  equispaced  samples  of  one  of  the  quadrature  components .  The  corre¬ 
sponding  sample  mean  is  subtracted  from  each  member  in  each  data  set, 
and  the  result  is  then  prewhitened.  Prewhitening  of  sample  data  im¬ 
proves  the  results  of  the  spectral  averaging  which  is  performed  later 
in  the  program. Let  W(t)  be  the  result  of  prewhitening  U(t);  then  the 
prev'hitening  characteristic  used  in  the  program  relates  W(t)  and  U(t) 
in  the  following  simple  manner: 


W(ti)  =  U(tj.)  -  0.6U(ti_].) 


(13) 


This  prewhitening  characteristic  multiplies  the  input  spectrum  by  the 
quantity3 


1.36  -  1.20  cos2nT\ 

where  T  is  the  sampling  period.  Observe  that  the  high  frequencies  are 
emphasized  relative  to  the  low  frequencies. 

The  next  step  in  the  power  spectrum  program  computes  the  sample 
autocovariances4;5  and  crosscovariances 


WV  '  Ruv(V 


w 


and  Ryu(Ti) 


where 


N-i 

(T  )  =  —  E  UCt  )U(t.  +  T.) 
i  N  j_i  j  j  i 
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(Ti) 


U(t.)V(t  .  +  T.) 
J  3  1 


and  similarly  for  R^(t^)  and  Ryu^Ti^  •  These  covariance  functions  are 
computed  for  lags  out  to  one- third  of  the  record  length.3 

The  results  of  the  covariance  computations  are  then  used  to  com¬ 
pute  the  one-sided  cosine  transform  of  l/2(Ryy(T.)  +  and  the 

one-sided  sine  transform  of  1/2  ^Ryy(T\  )  ~  Ryy(T\  )j  •  These  transforma¬ 
tions  are  consistent  with  Eqs.  (7)  and  (9). 

Spectral  averaging  is  then  performed  on  the  results  of  the  cosine 
and  sine  transformations.  Spectral  averaging  is  necessary  because  the 
Fourier  transform  of  the  sample  covariance  function  will  not  converge 
to  the  true  power  spectrum  even  for  an  infinite  number  of  samples.4 
This  computation  may  be  viewed  as  obtaining  an  estimate  of  the  average 
of  the  true  power  spectrum  in  a  small  region  about  a  particular  fre¬ 
quency  as  opposed  to  obtaining  a  point  estimate  of  the  true  power 
spectrum  at  that  frequency.  The  function  that  defines  the  small  region 
over  which  the  spectrum  is  averaged  is  called  a  spectral  window .3j4>5 
The  program  uses  the  hanning  spectral  window3  which  performs  the 
following  simple  averaging  on  the  results  of  the  cosine  and  sine  trans¬ 
formations: 

A(0)  =  1/2  P(0)  +  PC1/2NT) 

A(k/2NT)  =  ^  P(k  -  1/2NT)  +  |  P(k/2NT)  +  |  P(k  +  1/2NT) 

for  1  S  k  S  S  -  1 

A(1/2T)  =  1/2  P(N  -  1/2NT)  +  1/2  P(1/2T) 

where  P(k/2NT)  is  the  result  of  the  cosine  or  sine  transformation  at 
frequency  k/2NT,  and  A(k/2NT)  is  the  resulting  spectral  average  at 
frequency  k/2NT. 
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The  program  then  corrects  for  the  prewhitening  initially  performed 
on  the  input  data.  This  is  accomplished  by  merely  dividing  the  spectral 
average  A(k/2NT)  by  the  quantity 

Jr 

1.36  -  1.20costt— 

N 

Finally,  Doppler  profile  estimates  are  obtained  for  positive 
Doppler  frequencies  by  summing  the  averaged  and  corrected  sine  and  co¬ 
sine  transformations,  while  estimates  for  negative  Doppler  frequencies 
are  obtained  by  subtracting  the  averaged  and  corrected  sine  transfor¬ 
mation  from  the  averaged  and  corrected  cosine  transformation.  See 
Fig.  1  for  a  flow  chart  of  the  power  spectrum  program. 

The  power  spectrum  program  was  operated  in  two  modes,  a  1-minute 
mode  and  a  2-minute  mode.  The  1-minute  mode  processes  1,200  samples 
of  amplitude  and  phase  spaced  at  50-msec  intervals.  This  yields  a 
20-cps  bandwidth  with  a  power  spectrum  estimate  every  1/40  cps.  The 
2-minute  mode  processes  1,200  samples  of  amplitude  and  phase  spaced 
at  100-msec  intervals.  This  yields  a  10-cps  bandwidth  with  a  power 
spectrum  estimate  every  1/80  cps. 

Figure  2  displays  the  results  of  processing  a  test  case  in  which 
input  data  possessing  tones  at  1  and  3  cps  were  simulated.  The  test 
case  illustrates  the  accuracy  and  resolution  of  the  program.  Examples 
of  Doppler  profiles  measured  on  the  Fort  Monmouth  to  Palo  Alto  path 
(4100  km)  and  the  Thule  to  Palo  Alto  path  (5050  km)  appear  in  Figs .  3 
and  4.1  Amplitude  and  phase  measurements  were  obtained  on  both  paths 
by  transmitting  a  phase-stable  CW  tone  at  7.366  Me.1 

Since  the  purpose  of  this  memorandum  is  merely  to  explain  a  data 
processing  technique  used  to  compute  channel  Doppler  profiles,  the 
physical  interpretation  of  the  measured  Doppler  profiles  is  deferred 
to  a  final  report  to  be  issued  under  this  contract. 
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FIG.  2  DOPPLER  PROFILE,  1-MINUTE  MODE,  TEST  CASE 
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FIG.  3  DOPPLER  PROFILE,  1-MINUTE  MODE,  FORT  MONMOUTH, 
11  FEBRUARY  1964,  1424Z 
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FIG.  4  DOPPLER  PROFILE,  1-MINUTE  MODE,  THULE,  10  MARCH  1964,  0210Z 
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